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Computational-Experimental Processing of Boride/Carbide Composites by Reactive 

Infusion of Hf Alloy Melts into B4C 

Arturo Bronson and Vinod Kumar, The University of Texas at El Paso 

Project Summary 

This project primarily investigated computationally the infusion of liquid elements from 
the Hf family (i.e., Hf, Ti and Zr) into a capillary to understand the effects of surface tension and 
viscosity during processing of melts with a packed bed of B 4 C to form boride-carbide 
precipitates. Although the ultimate goal of the research endeavor is to enhance significantly the 
oxidation resistance of ultrahigh temperature ceramic composites (UHTCC), the processing of 
these reactive metals and non-metals must be improved dramatically to control the eventual 
development of the oxide scale. The research effort used computational thermodynamics 
and computational fluid dynamics (CFD) to predict phase equilibria and melt flow into 
pore channels, respectively, with strategic experimentation through the following integrated 
research objectives: 

Investigate the use of surface tensions of liquid Hf, Sm, Ti, Y and Zr to predict their infusion 
into a B 4 C packed bed by computational fluid dynamics. 

Investigate the extent of reactive infusion with strategic experimentation of liquid Al-Sm 
alloys containing Hf, Ti or Y. 

Investigate the phases formed as the liquid metal flows into the packed bed of B 4 C. 

Although the alloy melt infiltrates as a result of the surface tension, viscosity, the liquid/B4C 
reaction, and the diffusion of components, the surface tension of the melt serves as the key 

parameter affecting the depth of melt penetration of liquid elements (i.e., Hf, Sm, Ti, Y and 
Zr) into pores. 

The research applied computational fluid dynamics (CFD) to simulate the flow 
characteristics of the melt from the Hf family of elements, as well as Y and Sm, through pore 
channels. The volume of fluid (VOF) approach with the momentum equations incorporating 
surface tension and contact angle terms was used to track the interface of fluid penetration into a 
pore. A pore-network modeling incorporating surface tension terms through contact angle and 
pore throat/nodes reconstructions is also considered for predicting dispersion. 

A modified Semlak-Rhines model was developed for liquid metal penetration into a 
packed bed and considers the contribution of the forces incorporating surface tension, viscosity, 
gravity and end-drag. For liquid Ti, the surface tension dominated the capillary flow at small 
times (< 1 ps) with the viscous force having a synergistic effect at longer times (> 10 ' 4 5 s). At 
small time (i.e. , < 1 ps range), the penetration rate reaches a maximum in the m/s range at 
approximately 1 ps, but the rate decreases with increasing time at which the rate of penetration 
compares favorably with the Semlak-Rhines equation. Though the surface tension is sometimes 
considered as the driving force in capillary flow, the contribution of the viscosity to the flow may 
dominate the surface tension-viscosity dissipation as the flow occurs at longer times or achieves 
a steady-state condition. 

Strategic experimentation involved Al-Sm and Al-Sm-Me (Me=Hf, Y or Zr) alloy melts 
infusing into a B 4 C packed bed at temperatures > 1973K (1700°C). The fluid dynamics of the 
Al-Sm-Hf and Al-Sm alloy melts was investigated through varying spatial pore channels and the 
elemental components formed carbide and boride phases (e.g., AI 3 BC). 


Bronson-Kumar (University of Texas at El Paso) 


1 



1. Introduction 

Ceramic composites and intermetallics have become the core materials to meet the 
strength and extreme temperature needed for gas-turbine engines and 
hypersonic aerosurfaces [ 1 ]. The successful integration of an ultrahigh- 
temperature, ceramic composite (UHTCC) for these applications 
requires careful and optimal consideration of processing, mechanics and 
oxidation. The research focuses on predicting the fluid flow of liquid 
melts from the mainly the Hf family (i.e., Hf, Ti and Zr) into a packed 
bed of B4C. During infiltration, reactions between the alloy melt and 
B4C will form an HfB2-HfC/Hf-Y-Ti composite, which will eventually 
develop an HfC>2-Y203-Ti02 scale, as configured in Figure 1 . The oxide 
scale must be configured with precipitates embedded and aligned to 
minimize the stress/strain distribution, but the key to the configuration is 
the predictability and control of the fluid flow processing at ultrahigh 
temperatures. For example, the phase equilibria of the pyrochlore phase 
(or Y2Ti2C>7) within an Hf02-Y203-Ti02 scale may form as layers of 
Y2Ti2C>7, HfTi04, cubic-HfCL, or tetragonal-HFC>2 as the alloy melt 
reacts with B4C. In addition, accurate chemical behavior of the Hf 
systems at extreme temperatures is needed to adequately predict 
performance. The stable and metastable conditions may adjust the final 
microstructure of the actual oxide interphase protecting the matrix at 
ultrahigh temperatures. 

On a previous study, the computational-experimental approach 
was used to investigate the effects of thermal mismatch and growth 
stresses on the ZrB2/ZrC/Zr-Si [ 2 ] system. However, the operational 
boundaries must be pushed by using rutile instead of silica to control 
vaporization [ 3 ] and Hf components rather than Zr compounds to attain better ultrahigh 
temperature performance. In addition, to acquire the ultimate Hf protective-scale system, the 
processing of the Hf melts with B4C must be better controlled. With an expertise in 
computational fluid dynamics (Kumar) and in high temperature reactions (Bronson), the 
collaboration investigated liquid alloy infusing into a packed bed of B4C. However, to 
improve the characteristics of a boride/carbide system, the research studied the processing 
science of toward controlling non-reactive infusion into B4C. 



Figure 1 -- Oxide 
Scale developed 
from the ZrC> 2 - 
Y 203 -Ti 02 System 
to protect a 
boride/carbide/met 
al matrix 


2. Research Effort and Findings 

2.1. Consideration of Infusion Driven by Surface Tension (or Surface Energy): 

The dynamics of liquid metal wetting across surfaces or into porous media encompasses a 
wide-ranging of material processing (e.g., brazing and reactive infiltration), as reviewed by 
Delanny et. al [4], Eustathopoulos [5] Saiz and Tomsia [ 6 ] and Liquid-metal reactive infiltration 
offers a methodology for near-net processing of metal or ceramic composites by penetrating 
liquid metal into a ceramic preform (e.g., Al into mullite [7] or Zr into B 4 C [ 8 ]). Although 
originally elemental liquids of Hf, Ti and Zr were infused into B 4 C, as reported by Johnson, 
Nagelberg and Breval [ 8 ] for liquid Zr, more recently Maheswaraiah, Sandate and Bronson [9] 
studied the infusion of a liquid Zr-Si alloy with B 4 C to form ZrC and ZrB 2 precipitates embedded 
within a Zr-Si matrix. To optimize the processing, a better understanding of the infusion of 
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liquid metal into a ceramic preform is needed to predict and to design ultimately the metal matrix 
composites or ceramic matrix composites. 

The dynamics of wetting for liquid metals on surfaces differ from organic liquids because the 
interfacial energies of the former are greater than the latter, as noted by Saiz and Tomsia [10]. A 
consequence of the interfacial energies is the significant spreading rates of liquid metals on 
metals and nonmetals. For example, Saiz, Hwang, Suganuma and Tomsia [11] determined that 
the spreading rates of Sn-Ag solders at a temperature of 723 K (996°C) after contacting Fe-42 Ni 
p02< 10' 19 kPa. Also, Grigorenko, Poluyanskya, Eustathopoulos and Naidich [12] reported 
spreading velocities of 0.5 to 1 m/s for Ni and NiPd melts on a-SiC, as well as NiPdTi and 
CuSnTi melts on vitreous carbon, at 1473-1773 K (1200-1500°C) in a vacuum. The dynamics of 
wetting of liquid metals may incorporate the effects of surface tension, adsorption of active 
species, and interfacial reaction though the mechanistic explanations have evolved from the use 
of a force balance in the Young equation as explained by Lupis [13] to applying the molecular- 
kinetic theory such as Benhassine [14] for Cu-Ag melts on Ni surfaces and Swiler and Loehman 
[15] for reactive surfaces of the Al-0 system. A stochastic mathematical framework was 
investigated by Delgado and Kumar [16] to investigate the effects of heterogeneity in the pore- 
structure. Delgado and Kumar also investigated a multiscale approach to couple the micro (pore) 
- physics with the macroscale investigation [17]. 

To predict the fundamental concepts involving the fluid flow of liquid metal driven by 
surface tension into a porous structure of a packed bed, Semlak-Rhines [18] developed 
expressions assuming flow through capillary parallel tubes. Martins, Olson and Edwards [19] 
analyzed the Semlak-Rhines expressions and selected an infiltration rate parameter consisting of 
the surface tension-force to viscosity ratio independent of the liquid density for evaluating the 
effect of alloy composition on the flow into porous media. 

Although the dynamics of reactive wetting would fundamentally incorporate the mechanisms 
noted previously, and even the effects of a temperature spike from exothermic interfacial 
reactions, the fluid flow into the packed must be understood first to determine the differential 
between non-reactive and reactive wetting kinetics. Hence, the primary purpose of the present 
study was to determine the fluid flow driven by the surface tension-viscosity dissipation as the 
liquid metal penetrates into the pores of a ceramic preform. The present study focused on the 
surface tension-viscosity dissipation assuming an unreactive pore surface, and left the dynamics 
of reactive wetting of a ceramic substrate for future work. 


2.2 Theoretical Basis 

The starting equation is given by conservation of momentum equation in the axial direction 
for a capillary, as given by Semlak and Rhines [18], a modified form of the differential equation 
by Ligenza and Bemstein[20], assumes rapid wetting without any reaction between the liquid 
and solid and is given as follows: 

^(nR 2 [ph + pg (/ = 2 nRacos( 9 ) - (p h + PgO- ~ h )) 

7 1 7 (dh\ 2 (1) 

-7 rR 2 gh(p - p g ) - -7 tR 2 P [—) 

The rate of momentum accumulation on the left side of the equation equals the summation of the 
forces on the right side for the terms involving surface tension ( a ) with the contact angle (0), 
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viscosity (jli), gravity and the end-drag, respectively. For the simplified case, assuming gas 
density (p g ) and viscosity (p g ) are negligible, Eq. (1) becomes 


d I , dh\ „ N (dh\ 

p — [nR 2 h — ) = 2 nRo cos( 0 ) — 8 nph —) 
dt \ dt) \dtJ 


— nR 2 ghp — -nR 2 p 



( 2 ) 


When the fluid flow in the capillary has fully developed which can occur at sufficiently small 
times, the rate of rise is being slowed enough that the rate of change of momentum and the end- 
Rcrcos0t2ju and rate of penetration, dhdt=RacodiU.(0)8/ut as developed by Semlak and 
Rhines[18]. However, Eq. (2) does not result in explicit solution when t -* 0 when the mass 
entered in the capillary tube and its velocity is zero. The surface tension term is assumed 
constant and independent of the depth of penetration of the liquid in the tube. This implies that, 
for the right hand of the Eq. (2) to be constant, the liquid must accelerate at infinite acceleration. 

To address the singularity, we assume that the meniscus has already been formed as 

shown in the Figure 2. The height (h) of the meniscus in Eq. (2) can be separated into two terms, 

h and ho. The modified Eq. (2) becomes 

d / _ N dh\ , (dh\ 

— I ( pnR 2 h + pV L ) — 1 = 2 nRo cos(0) - 8 np(h + h 0 ) 

2 . 1 2 f dk \ 2 (3) 

-p(nR 2 h + V L )g--pnR 2 i K —) 

VLis the liquid volume bounded by the horizontal plane passing through the lowest point of the 
fluid surface and the plane containing the contact line, and h 0 is the maximum height of the 
liquid above the lowest point of the liquid surface for the contact angle 0. The initial conditions 
are as follows: 

= 0 ( 4 ) 

t =o 


, an 

hlt=0 = 0] dt 


2.3 Results of Applying Liquid Elemental Properties from the Ti Family, Y and Sm to 

Modified Semlak-Rhines (SR) Equation: 

For liquid Ti flowing into a capillary with a 10 pm diameter, the velocity was calculated by 
using the Semlak-Rhines (SR) equation (i.e., Eq. 2) and compared with the modified-SR 
equation incorporating the inertia and end-drag terms, as shown in Figure 3. The thermophysical 
properties were acquired from the recent review of Paradis et al. [21] and the viscosity value 
from Ishikawa et al.[22] who updated the measured data measured primarily with electrostatic 
levitation, a containerless experimental method. The melting temperatures were the same as 
reported by Massalski et al. [23] except for yttrium. Paradis, Ishikawa and Koike [24] reported a 
yttrium melting of 1796 K (1523°C) as compared with the yttrium phase diagrams (e.g., Ti-Y) 
depicting a melting temperature of 1795 K (1522°C) though the temperature differential is within 
the error of the measured physical properties. The viscosity data was extrapolated to 2573 K 
(2300°C) for comparison of the velocity of liquid moving into the capillary at an isothermal 
temperature of the elements. 

The velocity was calculated at the melting temperature of Ti. For times greater than 10 ps, 
both equations give similar velocities. However, for times less than 10 ps, the SR velocity of the 
liquid infusion into the pore increases significantly and unrealistically, but the modified-SR 
velocity plateaus at 7 m/s and then decreases with decreasing time (as seen in the Figure 3). In 
the modified-SR equation, it is observed that all surface tension force is being used to accelerate 
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the mass (inertia) in the capillary as the time approaches zero. The relative importance of each 
term in Eq. (3) with decreasing time can be more easily compared by dividing each term by the 
surface tension term (i.e., 2jtRacos(0)). At t=10' 7 s, the inertia term (or the left side of Eq. (3) 
indicating the summation of the forces on the right side becomes driven by the surface tension, 
as shown in Figure 4. The inertia term was not present in the standard SR equation resulting in a 
singularity. 

The surface tension force contributes towards accelerating the fluid, and quickly reduces to 
zero at around 10 ps as the end-drag (which is proportional velocity squared) quickly becomes 
significant. The end-drag term reaches a maximum at around the same time and accounts for 
Fdrag/Fa = ~0.9). The maximum velocity is also achieved at this time. The viscous term 
(which is proportion to the velocity as well as the depth of penetration) slowly becomes more 
important. Since the drag and viscous forces act as resistance and as more liquid enter the 
capillary, the viscous force increases, the liquid starts slowing down reducing in the drag force 
faster than the viscous terms. This results in slight increase in the inertia term peaking at around 
50 ps. After this, the viscous term becomes significantly more important than the drag or the 
FviscousFa=~l»Fdrag or inertiaFa ) resulting an approximation similar to SR (as 
supported by Figure 4). 


2.4 Discussion of Results Acquired from Modified Semlak-Rhines (SR) Equation 

The effect of the contact angle and surface tension on the maximum calculated velocities 
was evaluated further to determine the major parameter driving the infiltration of the elemental 
liquid metal into the pore. The contact angle of the liquid Ti meniscus on the inner wall of the 
capillary also has a significant effect on the velocity, as shown in Figure 5. The maximum 
velocity of the liquid Ti at its melting temperature decreases from 7 to 2.7 m/s for a contact angle 
increasing from 0° to 80°, respectively. Though liquid Ti penetration initially approaches greater 
than 1 m/s, the velocity of infusion appears to decrease asymptotically to less than 0.1 m/s with 
increasing times greater than a ms. However, Eq. (3) does not account for variation in the 
contact angle, which can change significantly as the liquid metal reacts with the surface and as 
Eoehman, Ewsuk and Tomsia [7] and Saiz et al. [10] have reported. The transient effect of the 
contact angle on the velocity of penetration is beyond the scope of this paper though the reactive 
wetting of the surface should consider the interfacial kinetics. 

The velocities of liquid Ti, Y, Hf, Zr and Sm were also determined with the modified-SR 
equation at their respective melting temperatures, as shown in Figure 6. For the group IV 
elements (Ti, Zr and Hf), their velocities especially the maximum velocity increases with the 
decreasing order of the atomic weight, but the trend does not follow with the inclusion of Y and 
Sm. The maximum velocity increases with increasing surface tension of Sm, Y and Ti, as shown 
in Table 1. However, the surface tensions for Ti and Zr have very similar values (Table 1) with 
only a Aa of 0.1 Nm , but the respective velocities differ with the differential, Av max , 
approximating 1 m/s. In addition, the Hf velocity with surface tension does not follow the 
increasing trend. The maximum velocity for Hf, Ti, Y and Zr were also calculated at a constant 
temperature of 2573 K (2300°C), but a similar behavior was observed without an increasing 
trend as shown in Table 2. 

The effect of the surface tension to viscosity ratio on the maximum calculated velocity 
was also evaluated, but does not have a trend at the melting temperature, as evident in Table 1. 
The synergistic effect of the surface tension and viscosity depends on the transient time, as 
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shown in Figure 5. However, in examining the terms of the forces (Figure 4), the ratio of the 
viscous to surface tension forces increases with increasing time. At an isothermal temperature of 
2573 K (2300°C), the surface tension to viscosity dissipation (or ratio, a/p) does indeed have a 
trend decreasing with increasing elemental melting temperature. Though the surface tension is 
sometimes considered as the driving force in capillary flow, the contribution of the viscosity to 
the flow may dominate the surface tension-viscosity dissipation as the flow occurs at longer 
times or achieves a steady-state condition. 


2.5 Conclusions from the Computational Analysis of Capillary Flow 

The surface tension-viscosity dissipation driving the fluid flow within a capillary, the high- 
temperature liquid metal infusion was analyzed for titanium, yttrium, hafnium and zirconium to 
predict their penetration into a packed bed. A model of the dissipation considered the 
momentum balance within the capillary to determine rate of infusion, which modified the 
Semlak-Rhines model developed for liquid metal penetration into a packed bed assumed as a 
bundle of tubes mimicking the porosity of a packed bed. The modified model considers the 
contribution of the forces incorporating surface tension, viscosity, gravity and end-drag. For 
liquid Ti, the surface tension dominated the capillary flow at small times (< 1 ps) with the 
viscous force having a synergistic effect at longer times (> 10' 4 ' 5 s). At small time (i.e. , < 1 ps 
range), the penetration rate reaches a maximum in the m/s range at approximately 1 ps, but the 
rate decreases with increasing time at which the rate of penetration compares favorably with the 
Semlak-Rhines equation. Though the surface tension is sometimes considered as the driving 
force in capillary flow, the contribution of the viscosity to the flow may dominate the surface 
tension-viscosity dissipation as the flow occurs at longer times or achieves a steady-state 
condition. 
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Table 1-Physical Properties at Melting Temperature of Liquid Elements 


Element 

Melting 

Temperature 

K(°C) 

Surface 
Tension 
(^mp) Nm 

§Viscosity 

(Pmp 10 3 ) kgm'V 1 

(Pa»s) 

Density 
(Pm P 10' 3 ) kgm 3 

^mp/pmp 

ms " 1 

Hafnium 

2504 (2231) 

1.612 

7.1 

11.82 

227 

Titanium 

1941 (1668) 

1.555 

3.3 

4.17 

471 

Samarium 

1347 (1074) 

0.430 

2.29 

7.02 

188 

Yttrium 

1795 (1522) 

0.804 

4.54 

4.15 

177 

Zirconium 

2128 (1855) 

1.464 

4.74 

6.21 

316 


*Thermophysical data were acquired from Paradis et al.[21]. 

§Viscosity values were extrapolated from experimental data of Ishikawa, Paradis, Okada and 
Watanabe [22], 


Table 2-Physical Properties of Liquid Elements at 2573 K (2300°C) 


Element 

Surface Tension 
(o) Nm 1 

§ Viscosity (ju) 10 3 
kgm'V 1 (Pa-s) 

Density (p) 10 " 3 
kgm ' 3 

o/p 10‘ 3 
ms ' 1 

Hafnium 

1.60 

6.69 

11.8 

239 

Titanium 

1.44 

2.16 

4.03 

667 

Samarium 

0.341 

N/A 

N/A 


Yttrium 

0.765 

0.491 

3.99 

1558 

Zirconium 

1.41 

3.00 

6.09 

483 


§Viscosity values were extrapolated from experimental data of Ishikawa, Paradis, Okada and 
Watanabe [22], 


Table 3: List of Symbols and Their Units 


Parameter Description 

Symbol 

Unit 

Capillary Radius 

R 

m 

Temperature 

T 

K(°C) 

Liquid Viscosity 

P 

Pa-s 

Surface Tension 

a 

N/m 

Liquid Density 

p 

Kgm' 3 

Contact Angle 

e 

Degrees 

Acceleration Due to Gravity 

z 

ms" 2 

Height of lowest point from the liquid curved surface 

h 

m 

Height of the liquid above the lowest point of the liquid surface 

K 

m 

Equilibrium Capillary Rise Height 

he 

m 

Non-Dimensional Viscosity 

P 

- 

Non-Dimensional Liquid Volume Above Lowest Point of Liquid 

Interface 

v L 

- 

Parameter (X = (R/h e )V L 

a 

- 
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t = 0 s 

A 


R (1 - sin (0)/ 

T 

:os (0) 



[X 

V 

/ 

/ 

A 


4 

i* 




t> 0 s 

\ 



-si 


Figure 2 — Capillary tube is lowered onto the surface of a liquid that wets the tube. (Left) 
Liquid interface arranges to equilibrium capillary shape when the bottom of the tube 
contacts the liquid. (Right) Liquid moves into tube under the influence of capillary forces 



1.E-07 1.E-06 1.E-05 1.E-04 1.E-03 1.E-02 

Time, s 


Figure 3 -- Calculated velocity determined by the Semlak-Rhines equation at the melting 
temperature of 1941 K (1668°C) for liquid Ti with and without the inertia and end-drag 

consideration. 
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Figure 4 -- Importance of the various terms in Eq. (1) with respect to surface tension (F a ). 



Time, s 

Figure 5 — The effect of the contact angle (theta) on the calculated velocity (v) using the 
modified-SR equation for liquid Ti at 1941 K (1668°C) flowing into a 10 pm diameter pore. 
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Figure 6 — Velocity of elemental liquid for Ti, Y, Hf, Zr and Sm calculated with the 
modified-SR equation at their respective melting temperature. 


2.6 The capillary infusion into a pore calculated by computational fluid dynamics 
simulations compared to the Semlak-Rhines Equation 

For an understanding of the infusion of liquid hafnium (HfL), which simultaneously reacts 
with the packed bed of B 4 C, the present study approaches the behavior from the effect of surface 
tension or surface energy. The surface related parameters should examine the thermodynamic 
approach and Navier-Stokes approach though the thermocapillary effect (or Marangoni forces) 
should also be incorporated, because the exothermic Hfi 7 B 4 C reaction increases the temperature 
sometimes by 200°C. For a synopsis presented here, the surface tension and energy is explained 
with the effects of the fluid dynamics on the infusion. 

The Young-Dupre equation, which 
considers the forces acting on the point of 
intersecting lines for the solid/gas (a sg ), 
liquid/solid (ai s ), and liquid/gas (ai g ) interfaces, 
acquires the surface tension as explained by 
Lupis [13] and Gallois and Lupis [25] by 
measuring the contact angle, as depicted in 
Figure 7 and the following equation: 

o sg = ais + aig cos 0 (1) 

For two fluids (i.e., liquid and gas or two 
liquids), the force of the interface created with 
the solid and liquid minimizes the free surface 
energy by decreasing the surface area, as shown Figure 7 — Sketch of contact angle resulting 

in Figure 7. from forces 

Infusion of Hf-Y-Ti melts through a 
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B4C packed bed heavily depends on surface and interfacial tensions of the liquid alloy with the 
gas/solid surfaces and solid surfaces, respectively, as the liquid alloy permeates the pores via 
capillary flow, reacts with B 4 C to form HfB 2 and HfC with any remaining Hf alloy melt exposed 
to a temperature gradient. Although the term of surface tension (a, work/area), or work required 
to change the interface by a unit area, as defined by Wagner [26] has been used, one should note 
that Hondros [27] differentiates between surface energies (e.g., mJ/m 2 ) and surface tensions (e.g., 
mN/m). Lupis [13] also uses surface tension (force/length) though simply multiplying and 
dividing a length unit relate the energy and force terms, which are henceforth distinguished by a A 
and a F , respectively. The effects of surface adsorption on the thermodynamic stability is 
developed from the Gibbs free energy, G, as reported by Wagner [26]. 

dG=-SdT+VdP+adA+2 Mi^ (2) 

where symbols represent the usual terms such as entropy (S), volume (V), surface tension (a) 
and chemical potential (p). The oxygen potential would affect the adsorption in a 
thermodynamic system and related according the following equation derived from Eq. (2): 

f _ If do A \ (3) 

° KTVain P oJ T4 , k 

TOtif) would have a significant effect on the surface energy (g a ), especially since oxygen 
dissolves significantly in pure Hf. 

The research effort simulated the infusion into a pore by using computational fluid 
dynamics with Fluent software and compared its results with the Semlak-Rhines infusion 
equation [18] to acquire the depth of infusion (d) related by a F as follows: 


d = 


Ra F tcos 9 


2 V 


0.5 


(5) 


They derived the equation for a liquid rising into straight capillary tubes of radius (R), with a 
contact angle (0), the viscosity (rf) of the liquid. The depth of infusion increases with decreasing 
contact angle, but it appears to have a minimal effect with contact angles less than 20°, as shown 
in 4. 
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Figure 8 -- Depth of infusion of liquid yttrium according to the Semlak-Rhines equation 
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The Fluent software was configured to have volume elements simulating liquid hafnium 
(i.e., 0.1 mm x 0.1mm with 49,284 elements), liquid yttrium, or liquid titanium above a solid 
matrix containing one pore, as shown in Figure 9. A pore size had a 1 0 pm diameter and 1mm 
deep with 244,442 elements. A gas phase of CO maintained at one atmosphere (100 kPa) was 
also configured beneath the bottom of the pore, though as a fluid CO or air are approximately 
equivalent. Above the liquid, the pressure of the system had a zero AP and negating the gravity 
parameters to simulate the flow. 



Figure 9 -- Yttrium/CO volume fraction interface in a single pore using fluent 


Selmak-Rhines vs CFD for Hf Liq 



Time (s) 

Figure 10 -- Selmak-Rhines vs CFD for Hafnium Penetration 
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The depth of infusion calculated according to CFD simulation is compared with the 
Semlak-Rhines equation for liquid Hf at contact angle 10°, as shown in Figure 10. The 
corresponding calculations for liquid Y and Ti at contact angles of 10° are shown in Figure 11 
and Figure 12 respectively. For the liquid Hf simulations from Fluent, the trends are consistently 
below the Semlak-Rhines calculation though the trends seem to be converging as the time is 
increasing. The reasons for difference in the depth of infusion may result from the modeling 
parameters used in the Fluent software. Also, the assumptions made within the Semlak-Rhines 
equation can take an account for the differences in depth of infusion. We are also currently 
working on performing the grid-independence study. However, the Semlak-Rhines calculation 
indicates that the depth of infusion increase dramatically with a contact angle of 10° which is 
expected. 


Yttrium Selmak-Rhines vs CFD 



Time(s) 

Figure 11 -- Selmak-Rhines vs CFD for Yttrium Penetration 

For the infusion of liquid Y with a contact angle of 10°, the CFD simulations seem to 
correlate with the Semlak-Rhines equation as shown in Figure 11 as time progresses though the 
pore. Again, the CFD simulations remain in the same range of penetrating depth into the pore. 
However, the Semlak-Rhines equation does predict as expected, the deepest penetration of the 
three liquid metals is hafnium due to the ratio of surface tension to viscosity. In addition, all the 
times were calculated for very brief moments of infusion and simulations are expected for longer 
times. A similar result was acquired for a Ti melt flowing into a 10pm pore capillary as shown 
in Figure 12. 

The rate of penetration is calculated from the Semlak-Rhines equation and is evaluated as 
a function of temperature after a specific time, 100 ps, as shown in Figure 13. Each line 
represents a different size diameter of the pore. The rate of penetration prediction through the 
derived equation of Selmak-Rhines equation increases as the temperature increases and is 
consistent with the varying pore size. 

Flow through pore-network was modeled by solving a system of linear equations 
obtained from the conductance law (mass-conservation principle) at each node for all connected 
networks using quasi-steady state flow approximations. Conductance of irregular cross-section 
pore-geometry in Stoke problem is approximated by variation calculus. The basic idea of this 
variation method is to guess a test velocity function, which satisfies the given no-slip boundaries 
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conditions. This test function consists of parameters, which are determined by minimizing a 
functional that represents the total potential energy of the flow field. The resulting velocity 
function is the approximation to the exact solution of the Stoke equations by the variational 
method. 


Titanium Selmak-Rhines vs CFD @10° 
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Figure 12 -- Selmak-Rhines vs CFD for Titanium Penetration 
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Figure 13 -- Effect of temperature on the calculated rate of penetration yttrium for 

constant pore size. 
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Samples 


2.7 Experimental Study on the Infusion of Liquid Al-Sm, Al-Sm-Me and Hf-Y-Ti 
alloys into B4C: The previous processing technique 
of infusing liquid Zr-Si alloys into B 4 C was changed 
to study Hf-Y-Ti melts and Al-Sm-Me (Me = Hf, Ti, 
or Zr) melts into a B 4 C packed bed. A graphite 
enclosure was used for the processing of the 
composite, as shown in Figure 14. The samples were 
placed at the top part of the graphite crucible packed 
with a B 4 C bed maintained in the hot zone of the 
furnace (i.e., 1750°C). If needed, the hot zone could 
have two layers of crucibles with each layer having 
four crucibles containing the liquid alloy above the 
B 4 C packed bed. The CO generated from the residual 
oxygen within the enclosure will serve as a carrier gas 
pulling oxygen from the hot zone to the cold zone to 
maintain the oxygen potential of the Al-C-0 system. 

The oxidation potential as determined by 
Maheshwaraiah, Sandate and Bronson [9] was 
controlled by liquid A1 from 600°C to 1000°C with 
less vaporization of Al. Aluminum which melts at 
661°C and rests at the bottom of Zr 02 -Y 203 crucible 
within which the graphite enclosure was located. 

Towards the 



Graphite Enclosure 


Solid Graphite rod 
Graphite Plug/Cylinder 
Alumnium/AI203 filler space 
Zirconia Crucible 

Alumina Tube 


Figure 14 -- Sketch of graphite 
enclosure used to experimentally 
infuse Hf-Y-Ti melts into B 4 C 


Al-C-O, 1000 K 


bottom was the zone 
where the Al melt 

was placed to seal the 
crucible thus 

controlling the 

oxidation potential 
within the graphite 
enclosure serving as 
closed 

thermodynamic 
system. For example, 
the oxidation 

potential of 10' 47 atm 
will be fixed at the 
triple point of 

A 1 /A 1 2 0 3 /A 1 4 C 04 , as 
shown in Figure 15. 
To minimize the 

oxidation of the 



Figure 15 - Stability diagram for a Al-C-O system 


graphite, deoxidized 

Ar gas was passed through the alumina tube and subsequently bubbled through dibutyl phthalate 
to isolate the graphite from the air. The flow rates for the argon gas were measured with a flow 
meter with an accuracy of ± lmL. 
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Temperatures were measured with two Pt-20 wt% Rh/Pt-40 wt% Rh thermocouples with 
the bead located on the inner sidewall of the MoSb furnace chamber though another similar 
thermocouple was used to determine the 
temperature profile of the alumina tube 
containing the graphite enclosure shown in 
Figure 14. 

A major problem in the experiments was 
the disintegration of the sample by moisture 
reacting with the yttrium carbide. Although we 
attempted to determine a method for the 
formation of YB 2 or even YB 4 , the yttrium 
composition must be decreased to minimize YC 
or YC 2 formation. Master Ti-Y-Hf alloy melts 
were first prepared with compositions as 
indicated with X marks shown in Figure 16. 

The compositions chosen for experimentation 
marked on the red line (Figure 16) lie on the 70 
at. % line for Ti content. Figure 16 ~ Compositions of master alloys 

For processing with Hf, Ti, Y and Zr shown on the Hf-Y-Ti phase diagram, 
alloys, the oxygen potential must be controlled 

as evident in the stability diagrams for the Ti-C-0 and Y-C-0 systems superimposed in Figure 
18. Yttrium and Ti oxidize to form Y 2 O 3 and TiO at approximately 10 ' 24 and 10 ' 19 atm, 
respectively. The 

thermodynamic calculations for 
the present study were 
determined with Facts age 
software and database (version 
6.4). The pseudo-isopiestic 
technique successfully reacted 
the liquid alloys of the Ti 
family (i.e., Ti, Hf, and Zr) and 
Al-Sm system without 
oxidation, as evident from the 
microstructures without the 
formation of oxides (Figures 21 
and 26) in the next section. 

Although the velocity 
of liquid Ti, Hf and Zr infusing 
into a capillary were calculated 
by the concepts of fluid 
mechanics according section 
2.4, above, the experimental findings determined that a critical parameter needing further study 
is the contact angle. The liquid alloy must infuse the porous morphology of the B 4 C packed bed 
shown in Figure 17. Liquid alloys of Hf-Ti-Y were placed on a B 4 C packed bed and even 
substituting Sm for Y as a component, but the liquid did not penetrate into the pores of the 
packed bed at 1973 K (1700°C) and sometimes the reaction reaching a temperature greater than 



Figure 17 -- SEM image of pore morphology of B4C 
packed bed 


y 
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2045 K (1772°C) as evident from Pt melting and placed adjacent to the graphite crucible 
containing the Ti-Y-Hf liquid on a B 4 C packed bed. The Hf-Ti-Y calculated phase diagram and 
the microstructure of Figures 16 and 21 respectively, indicated that a liquid was possible but 
none of the liquid alloys shown in the Ti-Hf-Y ternary system infiltrated into the B 4 C. 

The alloy system was changed to the Al-Sm-Me (with Me = Y, Ti, Hf and Zr) to decrease 
the melting temperature, as well as keep the similar characteristics of Y with its substitution of 
Sm. By lowering melting temperature of the alloy composition, we were aware that the contact 
angle of the liquid was also a critical parameter, as Halverson et al.[28] reported on A1 
infiltrating B4C. They determined that the contact angle of liquid A1 decreased from 148° to 28° 
with increasing temperature from 1173 K (900°C) to 1573 K (1300°C), respectively. In addition, 
the time at temperature affected the contact angle though it was more dramatic at 1173 K 
(900°C) for which the contact angle decreased from 148° upon immediate contact with B4C to 
22° after 1450 minutes. Upon considering the temperature differential above the melting 
temperature of the alloys and the reacting time, the Al-Sm-Me system was decided to study the 
reactive infusion into B4C. 


Y-Ti-C-O, 1973 K 


0 < Ti/(Y +Ti) < 1, = 10 atm P(total) isobar 


^kctSage 1 ' 



Figure 18 -- Superimposed stability diagram of the Y-C-O and Ti-C-O systems at 1973K. 


In addition to oxygen potential, Ellingham diagrams best depict the stability of the 
carbides and borides compared to B 4 C, as shown in Figures 19 and 20. The carbides of Sm, Y 
and Ti are clearly more stable than B 4 C at temperatures greater than approximately 1473 K 
(1200°C). However, AI 4 C 3 may form upon cooling from 1973 K (1700°C) as the intersection of 
the B 4 C and AI 4 C 3 occurs near 773 K (500°C). As for the borides, AIB 2 is the only boride more 
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stable than B 4 C except when the temperature is greater than approximately 1873 K (1600°C) at 
which B 4 C is more stable than SmB 4 . The borides and carbides for elements of the Hf family 
(Hf, Ti and Zr) are clearly more stable than B 4 C as similarly to the oxides, which were reported 
by Maheswaraiah, Sandate and Bronson [9]. A comparison of the stability of carbides and 
borides of the elements used in the system is shown in the Figure 19 and 20 


Delta G (Ellingham) Delta G (Ellingham) 
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2.8 Typical microstructures from experiments 

Typical microstructures obtained after using the pseudo-isopiestic technique to control 


Figure 21 -- Microstructure of Ti-Y-Hf Alloy Figure 22 -- Microstructure of Zr-Y-Ti/B 4 C 

oxygen potential for preparing the highly reactive elements (i.e., Ti-Y-Hf) and Zr-Y-Ti reacting 
with B 4 C are shown in Figure and 22, respectively. The microstructure of Hf-Ti-Y alloy 
showed no formation of its oxides and the primary phases are the high temperature (3-Hf phase, 
dissolved with Y and Ti in the expected body-centered cubic structure (P-phase), and the low 
temperature hexagonal closed packed (a-phase). The remaining matrix consists of primarily the 
Ti-Y-Hf phase. For a Zr-Ti-Y melt reacting with B4C, the precipitates of boride and carbide 
within a Zr-Ti-Y matrix appear in Figure 22. 



Figure 23 -- Ti-HfB 2 -C >2 phase diagram indicating the oxygen potentials for the oxidation of 

TiB 2 and HfB 2 phases 
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Oxide precipitates were not found as an indication that the oxygen potential was 
sufficiently low without the formation of Y 2 O 3 . For the Hf component, the oxygen potential 
must be less than 10’ atm must be achieved to ensure no formation of Hf 02 (Figure 23), as 
shown in the microstructure of Figure 24. Hafnium does precipitate along with AFSm 
surrounded by the Al-Sm-Hf alloy though the expected AF,Hf was curiously not distinguished in 
the microstructure. 



Al 3 Sm 

Undissolved Hf 
Al-Sm-Hf alloy 


Figure 24 — 75Al-20Sm-5Hf alloy (content indicated in wt. %). 


The reaction of a liquid Al-Sm alloy with B 4 C caused the formation of A1 borides and 
carboborides, as shown in the microstructure of Figure 25 complemented with its line scan in 
Figure 26. The phases formed between Al-Sm and B 4 C are identified in Figure 27. On the Al-B 
binary phase diagram, the formation of AIB 2 occurs at 1229 K (956°C) via the following 
peritectic reaction: 


AIB 12 + Liquid = AIB 2 (5) 

With the formation of AI3BC, the cooling the sample may have transformed the ternary structure 
of AI3BC3-AI2.1B51C8-AIL into AI3BC. 

The SEM image labeled as AFSm phase (or h-A^Sm) with the hexagonal structure 
results from the following catatectic reaction: 

fLAl 4 Sm = SAl 3 Sm + Liquid (5) 

The Al-Sm liquid forming from the solid pA^Sm phase with a tetragonal structure forms as it 
cools to precipitate the h-A^Sm, as explained by Zhou and Napolitano [29] in their study of the 
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Al-Sm binary system. They also reported that the aAlnSm 3 phase is metastable when the 
6Al 3 Sm precipitates. 

The formation of carboborides and Sm borides were detected by the XRD spectrums with 
the alloys reacting with the B 4 C. The possible reactions taking place in the system are shown in 
the following equations. 


aAlB 12 T L ^ A1 3 B 48 C 2 T A1B 2 

(6) 

A1 3 B 48 C 2 T L ^^ A1B 2 + A1 3 BC 

(7) 

4A1(1) + 3C(s) -»A1 4 C 3 

(8) 

Sm(l) + B 4 C —* SmB 4 + C 

(9) 

Sm(l) + B 4 C —► SmB 4 + SmC 2 

(10) 

3SmB 4 (s) —► 2SmB 6 (s) + Sm(g) 

(11) 


Halverson et al. the ternary Al-B-C at 1273 K. At lower temperatures, 1141 K, the 
formation of AI 3 BC and AIB 2 exists and at temperatures >114IK AIB 2 , because of its lower 
decomposition temperature of 1165K is replaced by P-AIB 12 , which is also referred to as 
AI 3 B 48 C 2 [30]. Though the AI-B 4 C study detailed particularly the Al-B-C ternary system, the 
microstructures clearly indicated the coexistence of the AIB 2 , AI 3 BC and B 4 C. However, in 
conferring with the Al-B system, AIB 2 was precipitated after cooling from 1248K as Halverson 
has reported. Hence at temperatures greater than 1248K, the AIB 2 phase exists. In addition, B 4 C 
phases are either slightly surrounded by the AI 3 BC phase suggesting seemingly peritectic type 
reaction or the reaction between liquid Al and B 4 C to form AI 3 BC. 

The samarium-boron phase diagram exists with a few borides SmB 4 , SmB 6 and SmB 66 for 
which the peritectics of SmB 4 and SmB 6 are estimated. SmB 2 which is not shown on the phase 
diagram was only reported to form at high pressures (6.5Gpa) [31]. SmB 6 has a congruent 
reaction with the liquid with a reported temperature of 2580°C. Similarly, the SmB 6 with on 
reaction with the liquid is further leads to the formation of aSmB 66 with a peritectic temperature 
of 2150°C. 

For the Al-Sm-Me reaction with B 4 C, one must consider the Sm gas formation as shown 
in reaction 11. The XRD spectrums show the presence of the samarium hexaboride (SmB 6 ), 
which may be the result of tetraboride to hexaboride conversion as described by Sturgeon and 
Eick [32], The formation of SmB 6 results in endothermic reaction (AH= 1247 kJ @1973K) and 
the Gibbs free energy indicates a favorable reaction (AG = -311 kJ). However, in the present 
study Sm vapor was not observed though possible. 

Samarium forms a brittle dicarbide (SmC 2 ) though very limited research is available for 
samarium and its use in alloys. Although a proper samarium-carbon binary system has been 
reported, but a review by Gschneider and Calderwood [33] lists more phases (aSm, pSm, ySm, 
Sni 3 C, Sm 2 C 3 , aSmC 2 , pSnKy in a Sm-C system. The reaction of samarium with boron carbide 
results in the formation of SmB 4 and C as shown in the reaction 9. The enthalpy of this reaction 
at 1973K is -430 kJ indicating the reaction is exothermic, but above 1808 K the reaction 
becomes unfavorable as the Gibbs free energy of the reaction is positive. The effect of Sm, 
which behaves similar to Y, on the reaction with B4C to form Sm carbides and borides requires 
further study. 
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Figure 25 -- Microstructure of a Al-Sm-Hf with B 4 C (80-15-5 wt. %) sample with a line 

scan 
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Figure 26 -- Line scan analysis of the microstructure Al-Sm-Hf 
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Figure 27 -- Al-Sm alloy Al-Sm (74-26 wt. %) on B4C acquired with a scanning 
electron microscope using backscattered electrons. 
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2.9. Objectives Achieved: The objectives met during the research effort: 

1. The infusion of liquid Hf, Ti and Y was simulated with a simple model of computational 
fluid dynamics by using Fluent software. 

2. The pseudo-isopiestic technique was designed for investigation at 1700°C for the Hf-Ti-Y 
alloy melt reacting with B 4 C though compositions limit the reactive infusion even though 
sometimes attaining or exceeding 2045 K (1772°C) because of the exothermic reaction. 

3. Hf-Y-Ti master alloys have been prepared for use in melt infusion into B 4 C packed bed. 
The key for alloys with yttrium is the minimization of the formation of Y2C and YC2. 

4. Al-Sm liquid alloys infiltrated into a B 4 C packed bed to form A1 borides and A1 
carboborides at 1700°C. 

5. A modified Semlak-Rhines equation was developed to simulate the effects of the capillary 
flow of the liquid as time decreases, and a maximum flow into the capillary is achieved at 
0.1 to 1 ps for liquid Hf, Sm, Ti, Y and Zr at their respective melting temperature. 

2.10. Future Goals : The research tasks for the next phase of the research are as follows: 

1. The Fluent software will include more pores for simulating the effect of adjacent pores 
connecting with each other. 

2. Expand the extended finite volume code to incorporate the boride and carbide precipitates. 

3. The experimental tasks at temperatures of 2000 to 2500°C could occur once an induction 
furnace is acquired to determine the effect of more reactive liquid alloys (i.e., Hf-Ti-Y- 
Re). Rare earth elements (e.g., Ce, Sm, Lu, Yb) may form carbide and borides to enhance 
the mechanical properties of HfC-HfB 2 -Hf composite. 
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3. Accomplishments/New Findings 

Accomplishments -- The reactive processing of boride/carbide composites with a 
controlled oxygen potential was published the International Journal of Applied Ceramic 
Technology. The methodology consists of a thermodynamic system equilibrating the Zr system 
at temperatures reaching near 2600°C with an A1 system at approximately 700°C. The 
temperature gradient within a pseudo-isopiestic technique established the oxygen potential near 
700°C. The processing technique is extended to form yttrium boride in a HfB 2 -HfC/Hf-Ti-Y 
matrix, which would incorporate a boride complex in the composite with possible strength and 
oxidizing improvements. The methodology of controlling the oxygen potential successfully 
melted master 70Ti-Hf-Y alloys, but their infusion into B4C packed bed has been problematic at 
1700°C. With the analysis of the fluid dynamics for liquid Y, the temperature and pore size 
affects the penetrating rate as shown in Figure 13. 

Relevance to Air Force Mission — The computational fluid dynamics of the infusion of 
Hf alloy melts could optimize the placement of precipitates to improve mechanical and oxidizing 
properties for HfB 2 /HfC/Hf-Y-Ti at 1700°C. The effects of surface energies, diffusion of 
components within the liquid phase, and interfacial reactions could be use to optimize the 
reactive infusion process. With the use computational fluid dynamics in liquid processing of 
boride/carbide systems at 1700°C, the processing simulations of ceramic composites could 
predict the formation of stable and metastable phases to more extreme ultrahigh temperatures 
upon coupling of the modeling with thermodynamics. 


4. Personnel Supported 

Sanjay Shantha-Kumar and Arturo Medina were primarily supported with AFOSR funds 
for their doctoral and master’s degrees, respectively. Sanjay Shantha-Kumar performed the 
strategic experimentation on the reactive infusion of Al-Sm, Al-Sm-Ti, Al-Sm-Hf, Al-Sm-Zr and 
Hf-Y-Ti melts into B4C. He continues as a research associate at the University of Texas at El 
Paso after acquiring his doctoral degree. Arturo Medina computed the infusion by using Fluent 
software with his background in fluid mechanics and heat transfer. Arturo served as an intern to 
work on fluid mechanics and/or heat transfer at the Kirkland Air Force Base for the summers of 
2014 and 2013. He joined the Federal Aviation Administration as a mechanical engineer for 
their technical staff after acquiring his master’s degree. 

5. Publications/Interactions 

In 2012, Maheswaraiah, Sandate and Bronson [9] published on the reactive processing 
the Zr alloy melts with B 4 C with a controlled-oxygen potential. The continuation of the 
processing technique was presented at the Materials Science & Technology 2013 meeting in 
Montreal, Quebec, Canada at the session on Innovative Processing and Synthesis of Ceramics, 
Glasses and Composites. The title of our presentation was “Rate of Reactive Infusion of Zr 
and Hf Alloy Melts into B4C with Controlled Oxygen Potential and by Surface Tension 
Driven Flow.” The computational analysis of the capillary infusion was recently submitted to 
Metallurgical and Materials Transactions with a manuscript entitle, “High Temperature Liquid 
Metal Infusion Considering Surface Tension-Viscosity Dissipation.” 

The research group is collaborating with Dr. Erica Corral of the University of Arizona by 
preparing samples using a plasma-spark-sintering furnace. David Pham, a previous graduate 
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student who finished his master’s degree with Professor Bronson, is pursuing his doctoral degree 
in Professor Corral’s research group. 

Professor Bronson visited Air Force Research Laboratory on March 27, 2014 to exchange 
technical views with Dr. Allan Katz and his colleagues on ultrahigh temperature ceramic 
composites. Dr. Bronson also presented a talk entitled, “Reactivity of Zr and Hf Boride- 
Carbide-Oxide Systems at Ultrahigh Temperatures.” 

6. Patent Disclosures/Awards 

No patent disclosures have been submitted, yet, though the university is reviewing a 
patent application for the processing of the reactive metals to control the oxygen potential. 
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8. Appendices 

In this appendix, steps involved in solving the modified SR equation are presented. We start with the 
Eq.(3) with the appropriate initial conditions. 

d / „ N dh\ , (dh\ 

— I ( pnR 2 h + pV L ) — I = 2nRo co s(0) - 8np(h + h 0 ) (3) 
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, , 1 (dh\ 

—p(jzR 2 h + V L )g --pnR 

dh 

h\ t=0 — 0; — — 0 (4) 

at f—o 

Dividing Eq. (3) by Jt R 2 yields 

2 

p — ({h + RV l } — ) = 2jtcrcos(0)/R - Sp(h + h 0 )R~ 2 — - pg(h + Rh L ) - \p(—\ 
dt\ dt) dt \dt) 

( 5 ) 

with V L =V L /(jrR 3 );h L =h L /R 


Eq. (5) may be written in the form 

p^-l{h + RV L }^) = pg(h e -h)~ 8 p(h + \ )R~ 2 ^ 1 p{ ^ 


where h c is the equilibrium capillary height given by 

pgR 

Dividing by the density p, we get 


-\{h + RV L }—\ = g(h e -h)- 


8 p(h + hf) dh j / dh 


dty "’dt) ~ e pR dt \dt) 

As small time, the rate of penetration is small and hence contributions from the h () can be ingored. 
We introduce the following non-dimensional capillary height and time 

h =h/h ;t = tJgTh {9) 


Then gh — {h + (R/h e )V L }— =gh e (l-h) 


ZpdgK r dh 


pR dr 


h—-\gh e \ — 


Dividing by gh e , we obtain 


±{{h + (R/h e )V L } d ^\ = (\-h)-^h d ^-\\ d ^ 


dr dr 


where jU = jU/(pR 2 ^fgh~) 


Writing v for dh I dt and using the fact that 

d 2 h d 
dr 2 dh 

we have 

(h +«)-i(iv 2 )+fv 2 =(l-h)-8phv 
dh 

with a = RVJh e 
Eq. (14) can be written as 


(13) 

(14) 

(15) 
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-i{* +af 1 f}-{l-h)-$/llij2flh +af 1 
ah 

with f(h) = \v 2 

Eq. (12) is subject to the initial condition 

/(0) = 0 

Special Case - fi = 0 

In this case Eq. (15) may be integrated explicitly, using Eq, (17), to obtain 


(h + a) 512 f = f (1 + a)\h + a) 512 -a 5/2 )-%\h + a) 1 ' 2 - 

This yields 


a 7 ' 2 ) 


(15) 

(16) 

(17) 

(18) 


/ = f(l + a)-f(/j + a)-fi/(h +af' 2 (19) 

with fi = |«r 5/2 (l + j a) (20) 

Note that, as h rises from zero,/increases from zero to a maximum and then falls back to zero 

again, at the maximum value of h , Since a is very small in comparison with unity, we have, 
approximately, 

L,H + t« (21) 

A better approximation is 

1 ^ 


1 + ^a 


^max = J + JOT- f« 5/2 T-V 

5 5 5 {/(! + «)} 


5/2 




cr 


5/2 


1 + f or 


( 22 ) 


(l + a) 5/2 

Noting that A max is larger than the equilibrium value of A max = 1, we see that the system is 
underdamped. We now compute the maximum value of f and hence the fluid velocity. Let h = h 0 
when f is maximized. The value of h 0 is obtained by setting the derivative of f to zero. 


C --! + !/SAA l +a') , ' 2 -0 

ah h 

Hence fi /(/? () + Cl) = /- (// () + Cl) 


(23) 

(24) 


Substituting Eq. (24) back into Eq. (19), with h set to h 0 we get 

Now, from Eq. (19), 

(h 0 + a) 712 = \ a 512 (1 + f a) 

Then /z 0 =((l + fa)|a 5/2 ) 2 -a 

General /a 

In this case Eq. (15) needs to be solved numerically, subject to the initial condition in Eq. (10). To 
do this, we divide the range of h into intervals \h kA ,h k \; k— 1,2, ... N, with h k = 0, h N —h . Let 

fk =f(K) (27) 


(25) 


(26) 
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( 28 ) 


For h k < h< h k+l 
Eq. (15) is approximated by 

+ af 2 f\- (1 - m + af 2 -8 j£(h+ af 2 j2f t -,h t s h < h M (29) 

ah 

Eq. (29) is solved subject to the initial condition 

f(h k ) = f k (30) 

Note that Eq. (29) may be written in abbreviated form as 

f\k +a f 2 f}-{\-B t h)(h +af 2 -,h, sh <h m (31) 

ah 

with B k = 1 + 8//^2 f k (32) 

Integrating Eq. (31) over the interval [/^,/^ +1 ] yields the result 

(L+«) 5/ 7* + , -(^+«) 5/2 / t -fa+vxL+«) 5/2 -H(L+< 2 

” 5 (1 + Bk a )(K +a ) 512 ~T^k(K + a ) in (33) 
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